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We propose a new mechanism of negative magnetoresistance in non-magnetic granular materials in 
which electron transport is dominated by hopping between two nearest-neighbor clusters. We study 
the dependence of magnetoresistance on temperature and separation between neighboring clusters. 
At a small separation we find a negative magnetoresistance at low temperatures and it changes over 
to a positive value as temperature increases. For a fixed temperature, magnetoresistance changes 
from negative to positive when the cluster separation increases. The change of magnetoresistance 
AR/R can be more than 80% at low temperatures. 

PACS numbers: 72.80.Tm, 75. 70. Pa, 81.05.Rm 



Magnetoresistance in various transport mechanisms 
has been the subject of many studies. Negative mag- 
netoresistance (MR) has been observed and explained 
in many systems in the past two decades Q. Negative 
MR in dirty metals is related to weak localization phe- 
nomena H. Negative MR in Mott's variable-range hop- 
ping(VRH) region is also possible @ though backscat- 
tcring is unimportant there, because electron states are 
highly localized and the probability of backscattering is 
exponentially small. In the VRH case, it was pointed 
out that the ensemble average should be on the 

logarithm of conductance so that interference [^fj4| be- 
tween forward tunneling paths is important. The replica 
treatment || shows that pairing of tunneling paths is im- 
portant to the average. A magnetic field can introduce 
phases to tunneling paths. Thus, the pairing is weaken, 
the localization length increases, and negative MR is re- 
sulted. Recently, there are reports that negative MR was 
observed in A1/A1 2 3 j| and gold p] granular materials 
in which the nearest neighbor hopping is the main trans- 
port mechanism. One of the interesting features is that 
negative MR occurs only near the percolation threshold 
HD- It i s un hkely that these experimental results can 
be explained either by the mechanism of band conduction 
or VRH conduction. In this letter, we propose a novel 
mechanism for negative MR in a granular material or 
a quantum dot array containing non-magnetic elements. 
We show that level shifting, magnetic confinement, and 
quantum interference in a phonon assisted hopping pro- 
cess can produce very rich physics. 

There are three energy scales which play important 
roles in a magneto-transport in granular metals and 
quantum-dot arrays. They are the thermal energy kT, 
the level spacing 5, and the Zeeman energy. When the 
quantum dot sizes are small enough such that the typical 
level spacing is comparable to the thermal energy, only a 



few states near the Fermi energy contribute to hopping 
conduction. Therefore, if one can manipulate the rela- 
tive level positions around the Fermi energy, one can then 
change the relative importance of these states to hopping 
conduction, which, in turn, may change magnetoresis- 
tance substantially. Depending on the relative direction 
between the electron magnetic moment and magnetic 
field, an external magnetic field can either increase or 
decrease the energy of an electron state through the Zee- 
man effect. This level shifting may enhance the phonon 
assisted hopping and, thus, lead to a decrease of the re- 
sistance, namely, a negative MR in the nearest-neighbor 
hopping conduction. It is interesting to note that this 
mechanism of negative MR is completely different from 
that in Mott's variable-range hopping conduction [^P]. 

In the dielectric regime of a granular metal, metallic 
clusters are embedded inside an insulating matrix 
Electrons reside on these metallic clusters, and they can 
move inside the material by hopping from one cluster 
onto another. Therefore, it is important to understand 
how an external magnetic field B affects the electron hop- 
ping between two nearest-neighbor clusters. Consider 
a system with two clusters in the xy-plane centered at 
(d/2, 0) and (—d/2, 0), respectively. Each cluster is mod- 
eled by a parabolic confinement potential (This potential 
has been used to explain many experimental results Q). 
The one-electron Hamiltonian is described by 



H = -^-(p- -A) 2 + V(x, y) + n B aB, 
Ira c 



(1) 



with 



V(x,y) 



\muj 2 r ((x- f) 2 + y 2 ) x>0 
\mujf{{x + |) 2 +y 2 ) x < 



where m is the effective electron mass, and uj r and lui are 
the simple harmonic oscillator frequencies parameteriz- 
ing the right and the left metallic clusters, respectively. 



f 



The last term in Hamiltonian ([j]) is the usual spin Zee- 
man interaction in which /is is the Bohr magneton and 
er = ±1 correspond to spin-up and spin-down states, re- 
spectively. The larger u) is, the smaller a metallic cluster 
will be. When the two clusters are separated far away 
from each other, each cluster can be described well |8| by 
eigenenergies and eigenfunctions of a simple harmonic 
oscillator which can be analytically solved with eigenval- 
ues -E nj > = (2n + \j\ + l)h^Jujf r + uj 2 - jhuj c + fi B vB, 

where n = 0,1,2,..., j — 0, ±1,±2, a — ±1 and 
u> c = Be/(2mc). In the absence of a magnetic field, 
u c = 0, E n j is 2(2n + \j\ + l)-fold degenerate. The 
degeneracies are broken due to the Zeeman effect. 

In the hopping regime, thermal energy kT is much 
smaller than the Fermi energy fi, and the two clusters 
are only weakly coupled. Hopping can be regarded as 
an electron jumping from a state in the left cluster 
to a state ip2 in the right cluster. In the tight-binding 
approximation, the tunneling matrix element is [pi 



R cx 



f {A-x)ip* x i) 2 ] U=o dy, 



(2) 



where 4>o — is the flux quanta. In a uniform magnetic 
field B, with symmetric gauge A = (— By/2, Bx/2, 0), 
the spatial parts of tpi and ip2 can be expressed in terms 
of eigenstate ?/Vj °f a simple harmonic oscillator as |^] 



ip2 = exp(-i—A ■ r)ipn,j(r2,0 2 ), 



(3) 



(4) 



where Aq — Bdy/A, n and j are the quantum numbers 
of a simple harmonic oscillators, r\ and B\ are the po- 
lar coordinates of (rcos6> — d/2, r sin 9), and r 2 and 6 2 
are the polar coordinates of (rcos9 + d/2, r sin 6). The 
phase factors in Eqs. and are due to the mag- 
netic field. They will give the usual interference on the 
tunneling matrix element. We will see later that this in- 
terference together with level shifting and the usual mag- 
netic confinement can induce both positive and negative 
magnetoresistance. 

In a process that an electron loses completely its phase 
coherence after it hops from one cluster onto another (in- 
coherent hopping process), we can describe a granular 
system as a random resistor-network jiofl in which each 
pair of metallic clusters is replaced by a resistor. If the 
electric field is small such that changes of the Fermi en- 
ergy and eigenenergies are negligible compared with kT, 
then resistance between two clusters can be calculated by 



l,r 



exp( 



2kT 



(| El - E r 



Ei -n 



M I)), (5) 

where /i is the Fermi energy. Ei and E r are the energies 
of states ipi and tp r , respectively, including the contribu- 
tion from electron-electron interactions. The sum is over 
all states in both left and right clusters. Equation (|5|) 
takes into account of a process that an electron in state 
Ei (E r ) absorbs (emits) an phonon of energy E\ — E r , 
and hops to unoccupied state E r (Ei). The Fermi energy 
is related to the number of electrons N in a cluster. At 
zero temperature, /i ~ \/~NTiLO, where u> is the harmonic 
oscillator frequency for the cluster. In the calculations 
shown below, we consider a process that one electron 
hops from one neutral cluster onto another neutral clus- 
ter. Electron-electron interactions are taken into account 
through the capacitor effect, namely, an extra charging 
energy, e 2 /(nr), will be added on every eigenenergy of 
single electron Hamiltonian ([!]) when an electron is put 
on a neutral cluster, r is the size of the cluster, k is 
the effective dielectric constant. The Coulomb charging 
energy on a cluster of nano-meter can be much larger 
than the thermal energy. However, for an electron tun- 
neling between two neighboring clusters of similar sizes, 
the charging energy is close in value before and after tun- 
neling process, thus, the Coulomb blockade effect is not 
significant. This is very much different from a Coulomb 
blockade system where electrons tunnel from a lead to a 
quantum dot . 

Given uji, uj r , the effective mass m, and the magnetic 
field B, it is straight forward to compute R(B) numer- 
ically using Eq.@. To obtain some realistic numbers, 
we use, throughout this study, m = l.lm e , the effective 
electron mass of a gold metal. We assume the atom- 
atom distance to be 1.6 A which is suitable for a gold 
metal. Let us consider a system with the Fermi energy 
fi = ll.75fi.wo- This corresponds to about 144 electrons 
on a cluster. The diameter of the cluster is then about 
D = 38.4 A. We choose ujq in such a way that Hloq is 
equal to the typical level spacing of this kind of clusters, 
i.e. hu) Q = h 2 /(mD 2 ) = 4.7 meV. In order to eliminate 
possible artificial effects such as resonance tunneling, we 
assume that the two clusters are not identical by choosing 
lji = 4.7 meV/h and uj r = 4.5 meV/fi. The open symbols 
in Figure ^ are numerical results of the dimensionless MR 
AR/R = (R(B) - R(Q))/R(0) vs B when the cluster- 
cluster separation is about d = 80 A. At temperature 
T = 6 K, the MR is negative. It starts to saturate around 
magnetic field B = 6 T. The change in resistance is about 
70%, comparable to giant magnetoresistance (GMR) ob- 
tained in magnetic layer systems [ p^l3] |. At T = 12 K, 
MR is still negative, but the resistance change is substan- 
tially reduced to about 10%. At T = 30 K, the MR be- 
comes positive, and the resistance increases by about 10% 
at magnetic field around 7 T. These results can be un- 
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derstood from the following arguments. At very low tem- 
peratures (KT <§; huii, hw r ), only those states near the 
Fermi level, which lies in general between two discreted 
levels of a cluster, are important to the electron trans- 
port. In the presence of a magnetic field, at least two of 
these states approach to the Fermi level (one from below 
and one from above) , due to the breaking of time- reversal 
symmetry through the orbital Zccman effect. According 
to Eq.(||), resistance will be reduced because of hopping 
between the two states. With increase of temperature, 
more states will play significant roles in hopping con- 
duction. The energy shifting toward the Fermi energy 
becomes less important, and the phase interference on 
the tunneling matrix element (Eq.(||)) dominates. This 
leads eventually to a positive MR at high temperature 
when the magnetic field is strong. 
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FIG. 1. AR/R vs B at different values of temperature T. 
Open symbols are the data for a two-cluster system with 
H = 11.757kj , uji = uj = 4.7 meV/h, uj r = 0.965tJo = 4.5 
meV/7i, and d = 80 A. (circle) T = 6 K; (square) T = 12 
K; (up triangle) T = 30 K. Filled symbols are the data for 
a 10 x 10 cluster array. The harmonic oscillator frequency of 
each cluster is random with Gaussian distribution function of 
mean loo and the mean square deviation O.lwo- A* = 11-7571^0, 
loo = 4.7 meV/h, and d = 80 A. (circle) T = 6 K; (square) 
T = 12 K; (up triangle) T = 18 K. 

The behavior of a granular material is very similar 
to that of a two-cluster system. We also carried out 
a similar calculations on a 10 x 10 array system with 
the Fermi energy fi = 11,75Tiwq. Each cluster is still 
modeled by a harmonic potential with a random u. lu 
has a gaussian distribution with the mean value loq and 
the mean square deviation O.lwo- The separation be- 
tween two nearest-neighbor clusters is 80 A. The results 
of AR/R = {R{B) - R(0))/R{0) vs B for various tem- 
peratures are shown in Figure [l] by using filled symbols. 
The behavior is very similar to that of the case with two 



clusters. 

According to the phase factors in Eqs. (0) and (Q), 
one expects that interference on tunneling matrix ele- 
ment t\2 should increase with cluster-cluster separation 
d. We test this by studying a two-cluster system of the 
Fermi energy [i = 11.757iwo, u>i = loq = 4.7 msV/h, and 
u> r = 0.965wo = 4.5 meV/Ti. The dimensionless MR 
AR/R = (R(B) - R(0))/R(0) vs B with various d at 
temperature T = 6 K is calculated. The numerical re- 
sults are shown by the open symbols in Figure |2|. It shows 
clearly that MR changes from negative sign to positive 
sign as the separation increases from d = 80 A to d = 160 

A. 
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FIG. 2. AR/R vs B at T = 6 K for different values of 
cluster-cluster separation d. Open symbols are the data for 
a two-cluster system with n = 11.757kc>o, = Wo = 4.7 
meV/h, and uj r — 0.965o;o = 4.5 meV/h. (circle) d — 80 
A; (square) d = 107 A; (up triangle) d = 160 A. Filled sym- 
bols are the data for a two-cluster system with fj, — 11.35ftwo, 
uji = uio — 4.7 meV/h, and ui r = 0.965a>o = 4.5 meV/h. (cir- 
cle) d = 80 A; (square) d = 107 A; (up triangle) d = 160 A; 
(diamond) d = 187 A. 

The general behavior of MR should be robust to a par- 
ticular choice of the Fermi level. To demonstrate this, 
we calculate R(B) for a different value of the Fermi en- 
ergy ji = 11.35fiuo, but with the same u>i and u r . Filled 
symbols in Figure § are AR/R = (R(B) - R{0))/R{0) 
vs B at T = 6K with cluster-cluster separation d = 
80, 107, 160, 187 A. Comparing the results with those 
of the Fermi energy fj, = 11.757k^0i we see indeed that 
the qualitative features remain the same. MR undergoes 
a sign transition as the cluster separation increases from 
d = 80 A to d = 187 A. 

It is worth to stress that the nearest neighbor hop- 
ping conduction in homogeneous materials leads usually 
to a positive magnetoresistance Jl(J due to the destruc- 
tive quantum interference and magnetic confinement. 
However, we have used a simple model to demonstrate 
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that the level shifting can lead to a negative magnetore- 
sistance in the nearest-neighbor hopping conduction in 
granular materials. The negative MR was observed in 
recent experiments on AI/AI2O3 Q and gold granular 
materials . Our theory predicts a sign change of MR 
from negative to positive as temperature increases, and a 
similar sign change was also seen experimentally Q . Ac- 
cording to the theory, a sign change should also occur as 
cluster-cluster distance increases since then the quantum 
interference and magnetic confinement will play more im- 
portant roles in the hopping process. This is consistent 
with the experiments where the negative MR was only 
observed near the metal-insulator transition (near the 
percolation threshold) and on the insulator side MM . 

Our theory suggests a new type of non-magnetic ma- 
terials which can produce a GMR (> 80%). Unlike the 
usual magnetic systems in which the change of spin- 
related scattering in a magnetic field is responsible to the 
observed GMR [jl2|,[l3), the mechanism here is the level 
shifting due to the Zeeman interaction. There is no hys- 
teresis in the system since it is non-magnetic in nature. 
This property might be an advantage in real applications. 

We have used a symmetric confinement potential to 
model the quantum dots. The spherical symmetry gives 
rise to a degeneracy beside the spin degeneracy at zero 
magnetic field which will guarantee that there are always 
some states moving closer to the Fermi energy as mag- 
netic field increases. Closing of energy levels will enhance 
the phonon assistant hopping which, in turn, leads to a 
large negative MR. However, slight asymmetry will not 
destroy the above mechanism. In a system with asym- 
metric dots, the spin degeneracy survives, and the level 
closing will be gradual. Thus, the effects discussed in 
this work, such as GMR and sign change of MG will 
occur at a lower temperature. On the other hand, in 
a real granular metal, metallic clusters tend to form a 
spherical shape due to a large surface tension. There- 
fore, a metal with a large surface tension in an insulat- 
ing matrix is a good candidate to observe the transport 
property discussed here. Also, one expects that the nega- 
tive MR can be enhanced by annealing a sample because 
this process can make metallic clusters more symmetric. 
The nano-technology allows us to make artificial atoms, 
molecules, and quantum dot arrays via surface gates in 
heterostructures. It has been showed that quantum dots 
in such systems can be well described by symmetrical 
potentials ||[[4|. This may be an ideal system to test 
our theory with tunable dot size and dot-dot separation. 
Further studies along this direction should be interesting. 
Although we only studied a specific 2D model, but the 
physics discussed in this paper is expected to carry over 
to other 2D as well as 3D models. 

In summary, we propose a novel mechanism of GMR 
in non-magnetic granular materials in which the main 
transport process is hopping between nearest-neighbor 
clusters. We find that magnetoresistance changes sign 



as the temperature or separation between clusters varies. 
The unexpected negative MR is due to the level shifting 
of the Zeeman effect. 
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